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Abstract 

We investigate dynamical stability of the ground state against a time-periodic and 
spatially-inhomogeneous magnetic field for finite quantum XXZ spin chains. We use 
the survival probability as a measure of stability and demonstrate that it decays as 
P(t) oc t~ 1//2 under a certain condition. The dynamical properties should also be related 
to the level statistics of the XXZ spin chains with a constant spatially-inhomogeneous 
magnetic field. The level statistics depends on the anisotropy parameter and the field 
strength. We show how the survival probability depends on the anisotropy parameter, 
the strength and frequency of the field. 

1 Introduction 

Dynamical stability of quantum systems was studied from various viewpoints, e.g. from 
those of quantum chaos and energy diffusion in general [TJ [51 [3] . In the field of quantum in- 
formation, a lot of attention has been paid to the fidelity, which also describes dynamical sta- 
bility [IE] BE HIS]. Thc fidelity is defined as F(t) = \(xp(t)\ip 5 (t))\ 2 , where \ip(t)) = U (t)\ip) 
and \ips(t)) = Us(t)\if>) are unperturbed (S = 0) and perturbed (S =fi 0) time evolutions, re- 
spectively. When the initial state \ip) is an eigenstate of the unperturbed Hamiltonian, F(t) 
coincides with a survival probability. The fidelity and the survival probability are also stud- 
ied in the context of quantum irreversibility [HI [101 U • Those studies on fidelity, survival 
probability and energy diffusion have so far concentrated on highly-excited states or coher- 
ent states in random matrix models and kicked spin models. In real physics of condensed 
matter, however, it is more essential to elucidate the dynamical stability of the ground state 
for many-body quantum systems, which are non-integrable, against perturbations. 

The properties of dynamical stability can depend on the level statistics of systems. Re- 
cently, the level statistics of XXZ spin chains with a time- independent random magnetic field 
is studied in Ref. [TJ]. When the anisotropy parameter A is nonzero, Gaussian orthogonal 
ensemble (GOE) level statistics appears for a large field although Poisson-like level statistics 
is seen for a small field. When A = 0, however, the level statistics always shows Poissonian 
behavior. Therefore, it is desirable to know how a time-periodic random field will affect the 
dynamics of quantum spin systems with and without GOE level statistics. 

In this paper, we will study the dynamical stability of quantum XXZ spin chains against 
perturbations induced by a time-periodic and spatially-inhomogeneous magnetic field. We 
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use the survival probability as a measure of the dynamical stability and show how the 
survival probability is related to the level statistics. By theoretical and numerical analysis, 
we will demonstrate that normal diffusion occurs in the linear response region. Numerical 
analysis will also exhibit how the survival probability depends on the anisotropy parameter 
A, the field strength and the field frequency. This work will complement the latest study 
on the energy diffusion in a frustrated spin system with a non-random magnetic field [13] . 

The organization of the paper is as follows: In Sec.[2j we introduce the model and briefly 
describe our numerical method. In Sec. [3l the survival probability P(t) is investigated 
in an analytical way. It is suggested that P(t) should show a power-law decay on some 
assumptions. Numerical results will be shown in Sec. |4] The parameter dependence of P(t) 
is discussed. Conclusions are given in Sec. [5] 



2 Model and method 

The Hamiltonian of the quantum XXZ spin chain on L sites under the time-periodic and 
spatially-inhomogeneous magnetic field is given by 

H(t) = H + H 1 (t), (1) 

where 

L 

Ho = J + SjSj +1 + ASjSj +1 ), (2) 

3=1 

L 

Hi (t) =J2 B jSm{Lut)S]. (3) 

Here, S" — (l/2)cr" and (crj, crj, cj) is the jth-site Pauli matrix; the periodic boundary 
conditions are imposed. The parameters J and A are nearest-neighbor exchange interactions 
and the anisotropy parameter, respectively. The local field strength Bj's are uncorrelated 
random numbers that obey a Gaussian distribution with the following average and variance: 

(B 3 ) = 0, (4) 
{B 3 B k ) = Bl5 3k . (5) 

Before calculating the time evolution, we note the symmetries of the model. In the 
absence of the perturbation \H\(t) = 0], cigenstates are classified by z component of total 
spin 5 t z ot (= J2j=i total wave number (K), parity (in the special cases of K — 0, tt) and 
spin reversal (for S£ ot = 0). Eigenstates with different symmetries are uncorrelated. When 
TL\{t) 0, only S* oi is conserved, and manifolds with different wave numbers, parity or spin 
reversal become mixed. Here we choose the S* ot = manifold. It should be noted that no 
transition can occur if the magnetic field is homogeneous. 

Let us define the survival probability as 

p(t) = (\{m\m)\ 2 ) avc = (p(t)) avo , (6) 

where () aV o means the average over the random magnetic field. We numerically calculate 
P(t) for 100 samples of spin dynamics that starts from the identical initial state but is 
perturbed by different sets of -B/s. The average of P(t)'s over those samples gives P(s). 
We choose the ground state of Ho as the initial state |^(0)}. \^(t)) is a solution of the 
time-dependent Schrodinger equation 

ih^-Mt)) =n(t)m)) = [Ho+mmm)}. (7) 
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In our numerical calculation, we take h = 1 for convenience. The solution consists of a 
sequence of the infinitesimal processes as 

\iP(t)) = U(t; t - At)U(t -At;t- 2At) . . . U(At; 0)|^(0)). (8) 

To calculate a time evolution operator U(t + At; t) for each time step At (= 10~ 3 ), we use 
the fourth-order decomposition formula for the exponential operator [131 114] . Our numerical 
calculation below is mainly carried out on the system of L = 12, whose S* ot = manifold 
involves 924 levels. Here we note that there is a rather large energy gap between the ground 
state and the first excited state because of the finite system size, while there is no gap in 
the unperturbed system for |A| < 1 in the thermodynamic limit. 



3 Theoretical analysis 

Now we analytically investigate the short-time behavior of the survival probability. While 
applying the method cultivated in the context of energy diffusion [T], we will choose the 
ground state as an initial state and take the average over the random magnetic field. Let us 
expand the wave function \4>(t)) in an adiabatic basis: 

3 

Here, 

1 



',•(*) - - / dt'E^t'), (10) 



rt 

and Ej (t) is the eigenvalue of the j th state at time t: 

H(t)\j(t))=E j (t)\j(t)). (11) 

Let us notice that Tt(t) = Ho at t = nT, where n is an integer and T = 2it/uj. The 
eigenvalues and eigenstates of the unperturbed Hamiltonian are supposed to be well solved: 

Wo|&>=eil&>. (12) 
Substituting Eq. into Eq. ([7]) and multiplying (k(t)\ on the left, we obtain 

Here we notice that the phases of the eigenstates are chosen so as to satisfy (j(t)\4r \j(t)) = 0. 
As Hi(t) = Zsmtot, 

a k {t) =uco^t ^ ^^ e^-'^it). (14) 

Integrating Eq. (fl4)) . we have 

a k (t) =a fe (0)+w f dscosujs V Z k3 (s)e^ ek ^~ 9 ^ s ^a :j (s), (15) 
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where Z k j(t) = (k(t)\Z\j(t)) /[E k (t) — Ej(t)]. Before deriving the survival probability, we 
need to calculate |afc(t)| 2 . 



\a k (t)\ 2 = |a fc (0)| 2 + 2Re 



Jo 



+ uj 2 I dsi / dsa coswsi coswsa ^ ^ Zk ] (si)Zl j ,{s 2 ) 

xe i ^ s ^- e ^- ek ^ + ^'^a j (s 1 )a*,(s 2 ). (16) 
Using Eq. IT5|) in the second term of the right hand side of Eq. (fT6|) , we obtain 



M0)| 2 + 2Re 



u I dsicosws! J2 ^(si)e 4[efc(si) -^ (si)1 |a fe (0)| 



+2Re 



I dsi / ds 2 coswsi cosws2 ^ ^ Zkj(si)Zjj>(s 2 ) 



xe *( s i)-^( S x)+^W-e i /( S2 )] % .,( S2)a * (0) 
+w 2 / dsi / ds 2 cos ujs\ COSUIS2 ^ ^ Z k j(si)Zl J ,{s2) 

Now let us define an occupation probability distribution: 

2 



P fe (0EE(|(^|^))| 2 ) avc = 



-iSj (t) 



W)) 



(17) 



(18) 



In particular, Pk(t n ) = (|afc(in)| 2 ) ave when i = t n = nT. The survival probability P(t) 
is defined by Pi(t) when the initial state is \<pi). In other words, the initial condition is 
Pfe(O) = Ski- From Eq. (j4|), we can say 



(Zkj{t))a,vc — 0. 



(19) 



Considering Eq. ([5]), we assume that 



(Z k] {t)Zi, 3 ,{t')) mc = B 2 6 kk ,6 jf f(k,j)g(t)5(t - t'). 



(20) 

Here, f(k,j) ~ Pk,j{£k ~ £ j)~ 2 j ancl PfcJ i s supposed to correspond to the transition proba- 
bility between \(f> k ) and \<f>j). The factor (e/. — £j)~ 2 comes from the denominator of Z k j(t). 
We considered that the time average of Ej(t) is about Sj. And g(t) is a periodic function 
whose frequency 2u since Z k j(t) should be a function with a frequency ui. Using Eqs. (|19p 
and ([2"U]). P k {t n ) is approximately given by the following for k = I: 

P k (t n )~l+cu 2 B 2 f" dscos 2 Lus /(fc>iMs){(|ai(s)| 2 >ave-<|afe(s)| 2 >ave}- (21) 
Jo it**) 

Here we assume that the integrands of Eq. (|2ip have almost constant values from t n to 
Then we obtain 



P k (t n+1 ) - P k (t n ) = uj 2 B 2 1 T ^ 



PkJ 

[£j — e k ) 



-\P 3 {t n ) - P k (t n % 



(22) 
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where 7 is a constant. Let us consider a continuous distribution P{x, t). Then we may apply 
the following correspondences: Pk(t n ) — > P(x,t) and Sj — £k — > Ax. If jJfc+i,fc/(£fc+i — £fc) 2 
is large enough, compared with pk+2,k/ {sk+2 — £k) 2 , Eq. ([22]) can be approximated by a 
diffusion equation 

gg(M) _ „ d 2 P(x,t) 

dt ~ dx 2 ' [6) 

where D — oj 2 B 2 ^. The solution of the diffusion equation is generally given by P(x, t) — 
{ATTDt)^ 1 / 2 exp(— x 2 /4Dt). Therefore, if the assumptions described above is valid, the sur- 
vival probability P(t) [= P(0, t)] should show a power low decay ~ i" 1 / 2 for t — t n = nT. 

However, such a normal diffusion is expected to occur only for the linear response regime, 
where the perturbation is not very large and the perturbation theory is valid. If the pertur- 
bation is very strong, Eq. ([23)1 cannot be derived from Eq. ((22)) because Pk+i,k/{£k+i — £k) 2 
for i > 2 is not negligible under large perturbation. In the non-perturbative regime, where 
the perturbation is so strong that the perturbation theory fails, P(t) is expected to decay 
faster than in the linear response regime. 

The spectral properties of the Hamiltonian should also be related to the survival prob- 
ability. When the level statistics shows GOE, the smooth level density due to a lot of level 
repulsion guarantees the continuum approximation used for obtaining Eq. (|23p . In the case 
of Poissonian level statistics, however, levels tend to make clusters and Pk+i,k/(zk+i — £k) 2 
for i > 2 is not always negligible. Therefore, Eq. (|23|) is not appropriate for that case. It is 
expected that the survival probability here shows a different behavior from that of a GOE 
case. 



4 Numerical results 

Figure 0] shows log-log plots of P(t) at each period of time. In Fig. 1(a), where A = 1 
and lu = 1.0, P(t) for Bq = 0.3 decays as P(t) cx t" 1 / 2 . The behavior is consistent with 
the theoretical result in the previous section. For B = 2.0, however, P(t) decays very fast 
and goes down to a saturation value in a short time. The saturation value is estimated to 
be about 0.002 if the occupation probability spreads homogeneously over all levels. When 
t = nT, the number of levels is 489 but not 924, which is the number of the Hilbert space 
of the system, since there are many degeneracies. In Fig. IDJb), where A = 1 and u = 0.5, 
P(t) for Bq = 0.3 shows slower decay than P(t) cx < -1 / 2 . This case (lj = 0.5 and Bq = 0.3) 
may corresponds to the near-adiabatic case. Namely, energy diffusion is slow since the 
perturbation is small and slow. On the other hand, P(t) for Bq = 2.0 decays faster than 
P(t) cx t^ 1 / 2 at first and goes down to the saturation value. Those properties of P(t) seen 
in Figs.QJa) andQJb) are consistent with the expectation mentioned in the previous section. 
Namely, the survival probability behaves as P(t) cx t~ x / 2 in the linear response regime, and 
P(t) decays faster than t~ x / 2 in the non-perturbative regime. 

The behavior of P(t) in Fig.[4jc), where A = and lo = 1.0, is very different from that 
of Figs. IHa) andBJb). For both Bq = 0.3 and Bq — 2.0, P(t) rapidly decays and seems to 
saturate to a higher value than the saturation value in Figs.H^a) and 1(b). For A = 0, the 
value is estimated to be about 0.02 if the occupation probability spreads homogeneously over 
the all levels. We may say that an anomalous diffusion occurs for A = 0. The anomalous 
diffusion is responsible for the properties of energy levels. Namely, when A = 0, the level 
statistics is Poissonian, although GOE level statistics is observed for A = 1 [12 . When the 
level statistics is Poissonian, which implies the appearance of level clustering, it is expected 
that the occupation probability distribution does not broaden very much but forms a wave 
packet in energy space and moves to higher levels. The wave packet reflects like a soliton 
at the highest levels and moves back to lower levels and reflects again at the lowest energy 
levels since the system size is finite. Such behavior was suggested in energy diffusion for a 
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time (n=t/T) 

Figure 1: Survival probability P(t) for (a) A = 1, uj = 1.0, (b) A = 1, u = 0.5, (c) A = 0, 
u> = 1.0. The crosses (x) are for B = 0.3 and the solid circles (•) are for B = 2.0. The 
dashed line is just an eye guide for P{t) oc i -1 / 2 . 
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near-integrable system in Ref. [13j . Since we take the average over the random field such 
soliton-like behavior is also averaged. The behavior of the soliton-like wave packet is different 
in each trial. Therefore, the averaged survival probability has very rough data for A = 0. 

5 Conclusions 

We have investigated the survival probability P(t) of XXZ spin chains under spatially- 
random and time-periodic field. The survival probability decays as P(t) oc t -1 / 2 in the 
linear response regime. The property was derived theoretically and confirmed by numerical 
calculation. Numerical results also demonstrated that Pit) decays more slowly for small B 
or uj than i -1 / 2 and faster for large Bq or lu. When A = 0, however, P(t) decays and soon 
goes down to a saturation value. The difference between the behavior for A = 1 and for 
A = is caused by the property of energy levels, i.e. level-clustering or lcvcl-rcpulsion, for 
each A. 
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